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OS ' Abstract 
OS ' 

J>> A new version of differential renormalization is presented. It is based on 

pulling out certain differential operators and introducing a logarithmic 
dependence into diagrams. It can be defined either in coordinate or mo- 
mentum space, the latter being more flexible for treating tadpoles and 
diagrams where insertion of counterterms generates tadpoles. Within 
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this version, gauge invariance is automatically preserved to all orders 
in Abelian case. Since differential renormalization is a strictly four- 
dimensional renormalization scheme it looks preferable for application 
in each situation when dimensional renormalization meets difficulties, 
O ' especially, in theories with chiral and super symmetries. The calcula- 

tion of the ABJ triangle anomaly is given as an example to demonstrate 
simplicity of calculations within the presented version of differential 
renormaliz at ion . 
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1 Introduction 



Various versions of differential renormalization [1-7] explicitly describe how a prod- 
uct of propagators corresponding to a given graph and considered in coordinate space 
can be defined as a distribution on the whole space of test functions, starting from 
the subspace of test functions which vanish in a vicinity of points where the co- 
ordinates coincide. Differential renormalization provides a strictly four- dimensional 
renormalization scheme useful for applications f|, [|. It is based on 'pulling out' 
some differential operators from the initial unrenormalized diagram; these are either 
□ = d a d a [1-4,7], or S = \-^—x a || ||. The latter version looks more natural for 
generalization to the case with non-zero masses when one uses the operator 

S = - ^ t: — Us ^ 771;— , (1) 

2Y9ui ' 2 Y dnn' V 1 

where m are independent coordinate differences. 

The purpose of this paper is to present a new version of differential renormalization 
which automatically preserves gauge symmetry at least in the Abelian case. To do 
this we shall modify, in the next section, the version of ref. |5|, |J based on operator ([!]). 
The renormalization is defined recursively, for every graph, provided it was defined for 
graphs and reduced graphs with a smaller number of loops. It reduces to pulling out 
the operator S and inserting some logarithm into the given Feynman amplitude. The 
main idea of the new version of differential renormalization is to 'spoil' the diagram by 
the logarithms in a minimal way, in the sense that a minimal number of propagators 
is spoiled by the logarithms. 

In Section 3 we shall present a gauge-invariant differential renormalization scheme 
for QED. The underlying idea is simple: to insert the logarithms into the photon 
lines. This enables us to perform, for renormalized Feynman integrals, the same 
manipulations for establishing the gauge invariance as for unrenormalized ones. Thus 
gauge invariance is maintained automatically, to all orders.^ One-loop polarization 
of vacuum requires special treatment when an operator with explicitly transverse 
structure is pulled out before renormalization. 

In Section 4, by translating renormalization prescription into the language of mo- 
mentum space, we shall present general recipes applicable to arbitrary diagrams. In 
particular, we shall treat, in a more simple way (as compared to ref. |J), renormaliza- 
tion of tadpoles and diagrams where insertion of counterterms for subgraphs generates 
tadpoles (e.g. the sun-set diagram). Furthermore, the calculation of the triangle ABJ 
anomaly is presented as an example of the application of differential renormalization. 

Finally, Section 5 contains discussion of the obtained results. 

3 In Jl0| a recipe of differential renormalization formulated in momentum space and applicable 
for diagrams with simple topology of divergences was applied to prove relations relevant to Ward 
identities and corresponding to some partial classes of diagrams that contribute to Green functions. 
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2 Differential renormalization 
in coordinate space 



2.1 Renormalization prescriptions 

A Feynman amplitude is defined in coordinate space through products of propagators 

n r (xi, . . . , x N ) = n r (x) = Yl jGi(x n+ Q) - x n _(fi), (2) 

i 

taken over lines of a given graph T. Here 7r±(l) are respectively beginning and the 
end of a line /, and 

Gi(x) = — Pi{-td/dx,m) = Pi(-id/dx,m)G(x) (3) 

47T v — x 2 + zO 

is a propagator, with P/ polynomial and K\ a modified Bessel function. 

The strategy of differential renormalization is an explicit realization of the exten- 
sion of functional (0), from a subspace of test functions to the whole space (see details 
in refs. [||, Instead of renormalization prescriptions of 0, |6|, let us now define 
renormalization of Feynman amplitudes by the following recursive procedure. Note 
that we need to renormalize not only Feynman amplitudes themselves but as well 
Feynman amplitudes nf' z (x) 'spoiled' in some way by additional logarithmical depen- 
dence. This notation implies that the j-th power of some logarithm is introduced in 
a certain way into some fixed /-th line: 

nf''(x) = \G\ 3 \x^ +{l) - x v _ {t) ) II \Gi(xk + {i<) - x w _ m ). (4) 

I'^l 

For example, one can substitute Gi(x) (corresponding to this very line) by 

Gp } (x) =ln J >VG,(x). (5) 

However this is not the simplest version suitable for applications. A better variant is 
based on multiplication by a logarithm in momentum space: 

G?\x) = j±y 4 J d 4 g (-iy h? (H(q)/f) G l (q), (6) 

where the function H can be chosen as — q 2 — iO or mf — q 2 — iO. The latter ver- 
sion generally looks much more simple from the point of view of the calculational 
simplicity. In all these situations, the logarithmically spoiled propagator satisfies 

DGjp{x) = jGt 1} (x) - (1 + ai)Gjp{x) (7) 

with 
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and ai degree of the polynomial Pi in ([!]) . 

Suppose that we know how to renormalize all the Feynman amplitudes Ui f ,l (x) 

and TlpL(x) corresponding to reduced graphs r/7 and proper subgraphs 7 C T of 
the given graph T. (These subgraphs and reduced graphs have smaller number of 
loops, h. Recall that the reduced graph r/7 is obtained from T by reducing each 
connectivity component of 7 to a point.) Here / is a fixed line of T, and j' is arbitrary 
integer. This means that we know all counterterms contributing to the .R-operation 
(i.e. renormalization at the diagrammatical level) which acts on the corresponding 
Feynman amplitudes: 

MI^= £ A( 7 i)...A( 7i )n^ 

= i?TI^ + A( 7 )n^. (9) 

where 7 stands either for a subgraph or reduced graph, A is the corresponding coun- 
terterm operation, and the sum is over all sets {71, ... ,7^} of disjoint divergent 1PI 
subgraphs of 7, with A(0) = 1. The operation B! is called incomplete i?-operation. 

In this and the following section we shall consider diagrams without tadpoles and 
such that insertion of counterterms for subgraphs does not generate tadpoles. Note 
that this is sufficient for QED. 

Let us now define renormalization of the Feynman amplitude IBV as: 

ry^ 1 = _I_ (p + 00/2) (1 - m"- 1 ) #nf +M - -7——— J2 RC^ v +h \ (10) 

3 + J + 1 7cr . jg 7 

Here R' is incomplete -R-operation given by @, and u = Ah — 2L + a is the degree 
of divergence (with L the number of lines and a total degree of the polynomials Pi 
in (HD). The sum in the second term of the right-hand side of (|10D runs over all 1PI 
proper subgraphs 7 of T that do not include the line /. The operator M r performs 
Taylor expansion of order r in masses and external momenta of the graph T. 
Furthermore, 

A(r)nf'' = ifflf''- J R'nf' i . (11) 

Finally, the operations C 7 are determined from equations 

(s + u/2) Ru r = Y / c 7 i?n r = £ ^ c 7 n r + c r n r , (12) 

7cr 7cr 

where the operation C 7 inserts a polynomial P 7 of degree ^(7) in masses of 7 and its 
external momenta into the reduced vertex of the graph r/7. Symbolically we write 

c 7 n r = n r/7 o? 7 (13) 

where o denotes the insertion operation. In the language of coordinate space, 

C r U r (x)=V r (d/dx i ) J] S^\ Xi -x i0 ), (14) 
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where io is a fixed vertex. The polynomial for the graph Vr is expressed from the 
difference of the left-hand side of (0) and the sum in the right-hand side. 



2.2 Comments and examples 

1. A proof of the fact that relations (|10| ) -(|T2| ) define a correct /^-operation can be 
obtained by a trivial modification of the arguments presented in ref. || . In particular, 
it is proved that in the coordinate space with deleted origin (i.e. when at least one of 
the coordinate difference variables Xi — Xj is non-zero) the incompletely renormalized 
Feynman amplitude R! Ilp can be written as the right-hand side of eq. ( |I~0D . Then this 
equation enables us to extend the initial functional to the whole space of test functions 



so that, by definition, the right-hand side of eq. ( |i0|) determines the renormalized 
quantity i?IIp''. Indeed, the first term there does not have divergences because all 
sub divergences are removed by R' . As to the overall divergence, it is removed by the 
product of the operators (S + uo /2^j (1 — M 1 ^ -1 ). The preliminary Taylor subtractions 
result in producing a polynomial in momenta and masses of degree u and thereby 
reduce the degree of divergence from u to zero. After commutation of the operator 

+ u/2^ with this polynomial, the resulting operator S acts as in logarithmically 
divergent case (see details in || ||) by removing the divergence by multiplication of 
the first degree monomials contained in operator (|1|). 

The other terms in the right-hand side of eq. (|T0|) are manifestly renormalized 
Feynman amplitudes corresponding to reduced subgraphs (with smaller numbers of 
loops). 

2. Examples, (i) If a diagram is primitively divergent (i.e. does not involve 
sub divergences), with degree of divergence u, then the above prescriptions give 

ifflr = (§ + u/2) (l - M"- 1 ) nj?', (15) 

where I is an arbitrary line of V. Instead of Ilp'' one can also use a linear combination 
Z)^II r ' z , with — 1 ( see e -g- Sect. 4.2 where the triangle anomaly is calculated). 

(ii) To spoil the minimal number of the propagators involved it is natural to 
choose, as the line I, a line that is chosen for renormalization of some maximal sub- 
graph of the given graph. Let V be logarithmically divergent and let it involve only 
one proper divergent subgraph 7, with a; (7) = 0. Then, using some I G 7, one has 



ifflr 



~£ (U^SIL^) . (16) 



A generalization of this formula for the case when all the divergent subgraphs of 
the diagram form a nest 71 C 72 C . . . C 7 r = T, with = 0, and the initial 

Feynman amplitude itself involves the j-th power of the logarithm, looks like 

ifflr = |^y5 (n^S (n rV> _ 2 . . . SU^ 1 )) , (17) 
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where I G 7 i. 

(Hi) Let T be logarithmically divergent and let it involve only two disjoint divergent 
subgraphs with u>(ji) = 0. Let U G ji, i = 1, 2. Then one can define 

ifflr = sn rV7lU72) fen^) (stl^) - r (c 72 n r ^) , (is) 



.2 

where the operation C 72 is defined by 

C 72 n 72 = c 72 5^ ^(xj %i )- 

(iv) Let r possess the following family of logarithmically divergent subgraphs: 
71, 72, 712 = 7i H 72 and T itself. Let I G 712- Then 

iffl r = SOTj. 1 ', # = 1 + A( 7l ) + A(7 2 ) + A( 7l2 ). (19) 



Here A( 7i ) = i2( 7i ) - R'^), A( 7l2 ) = fl( 7 i 2 ) - 1 are found from (|T§) and (H) for 
J = 1. 

(u) Let T be as in the previous example and 712 = / (which is not of course diver- 
gent). Let us keep in mind the 2- loop photon exchange diagram of QED contributing 
to the vacuum polarization. We have uj(T) = 2 and w( 7 i) = 1 for one-loop vertex 
subgraphs. Then 

MI r = (S + 1)(1 - M^R'U 1 /, (20) 
R! = l + A( 7l ) + A( 72 ), A( 7i ) = R( 7i ) - 1, (21) 
i?n 7! = (5 + l/2)(l-M°)n 7 ;', (22) 

= i(5 + l/2)(l-M°)I^'. (23) 

3. We use here the operator (§ + ui/2j (1 — M 1 ^ -1 ) instead of a differential oper- 
ator of order uj + 1 of the previous version |5], [(J . Although this mixture of coordinate 
and momentum space operators might seem strange we insist that this version of 
differential renormalization is more simple that previous ones and this preliminary 
momentum subtraction is natural. In Section 4 this version will be translated into 
momentum space which happens to be appropriate for treating arbitrary diagrams. 
To apply the operator (1 — M w_1 ) one can turn to momentum space, perform Taylor 
expansion and then come back to coordinate space. However this momentum subtrac- 
tion is as well easily described in coordinate space. In particular, in the massless case, 
we have the following prescription for the action of the functional (1 — M^~ l )Il(x) 
on a test function 6(x): 



(0 



A^ _1 )n(x), <j){x)) = J dxU(x)(l-M^- 1 )(f)(x), (24) 
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where M^ _1 performs Taylor expansion in x. 

4. It is also possible to 'logarithmically spoil' the diagrams in the language of 
the a-representation, by inserting factors like In- 7 fi 2h D(a) into the integrand of the 
a-representation (here D(a) is a standard homogeneous function). However this pos- 
sibility looks disadvantageous because of the lack of control of manipulations that are 
relevant for establishing desired symmetry. Note that momentum space modification 
(H) of the propagator is easily described in the a-representation as insertion of In n' 2 oti 
(with \jt! proportional to \i in @). 

5. For the products of propagators in coordinate space all the vertex are con- 
sidered as external. Feynman amplitudes are obtained from the products Hp by 
integrating over coordinates associated with internal vertices: 

F r {xi, ...,x n )=J dx n+1 . . . dx N U r (x 1 , . . . , x N ). (25) 

When writing down renormalization prescriptions flTD|)-(|T2"D for Feynman amplitudes 



(|25|) one obtains similar formulae where the operator S is now given by the sum over 
external coordinates. 

6. The coefficients polynomials P 7 play the role of contributions of simple poles 
to counterterms within dimensional renormalization. They are related with renormal- 
ization group functions by the same formulae as in the case of the previous version 
of differential renormalization — see l6fl. 



3 Gauge- invariant differential renormalization of 
QED 

As is well-known the Ward identities in QED, e.g. 

q a T a (p,q,p + q) = S(p + g) - E(p), (26) 

that connects the vertex and self-energy Green functions are proved using standard 
manipulations based on the following identity involving the electron propagator: 

J- [S(x 2 - Ol*S(Z - Xl )} = -[5(x 2 - - 5(f - x 1 )]S(x 2 - x x ), (27) 

which are based on the equation of motion for the free electron propagator S(x) 
[in ■ i 0)G(x), namely (m — i 0)S(x) = S(x). 

In fact one uses (see, e.g., |TT|) a natural one-to-one correspondence between 
diagrams that contribute to these Green functions (vertex diagrams are obtained 
by insertion a new triple vertex into lines of the electron path between external 
electron lines). The main problem in establishing the Ward identity is to prove that 
these manipulations are also valid for the renormalized Feynman amplitudes. Let 
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us now use the structure of our renormalization procedure which reduces to pulling 
out differential operators and spoiling the Feynman amplitudes by logarithms. Note 
that commutation relations of the differential operators with monomials in external 
momenta (i.e. derivatives in coordinates) are very simple: 

p r S ={S + r/2)(l - M r - r )P r , (28) 

where P r is such a monomial of degree r (the value r = 1 is relevant for (p6|) ) . 

Therefore, the problem is not to spoil identities fl27|) by the inserted logarithms. 
Since these identities are connected with electron lines a natural and simple solution of 
this problem is to introduce the logarithms only into photon lines. Then the proof of 
(|26|) is performed by induction (as the renormalization prescriptions themselves), with 
the use of (|10|), under assumption that the corresponding relation between vertex and 
self-energy diagrams is valid for all subgraphs and reduced graphs. A crucial point is 
that in the right hand side of (|10| ) one has reduced graphs with at least one photon 
line if the initial graph has a photon line, because summation in the right-hand side 
is over subgraphs 7 such that /G7. 

The Ward identities that connect iV-point and iV — 1-point Green functions, with 
N > 3 and N = 2, are analogously proved. Only diagrams without photon lines 
require special treatment. These are the one- loop polarization of vacuum, vertex 
and box diagram. However, the second one is zero, the third is convergent. Thus, 
to complete consideration of the case N = 2 it is sufficient to prove that the one- 
loop polarization of vacuum is transverse. A straightforward application of general 
prescription (10) leads to a non-transverse expression. Of course, it is possible to 



adjust finite counterterms and arrive at a gauge-invariant result. (Why is it bad to 
do this just for one diagram of the theory?) However, following the style of differen- 
tial renormalization, one can modify the corresponding prescriptions by pulling out 
differential operators in a manifestly gauge-invariant way. Using manipulations |7] 
based on standard properties of Bessel functions, the one-loop polarization of vacuum 
can be written for x 7^ as 



2 / im 



2\ 2 



n, v (x) = -e l I — 1 tr{(m + * 0)f(z) lfl (m-i 0)f(z)} 



1 4 
em 



12tt 



(d^-g^)h(z), (29) 



where 



f(z) = K 1 (z)/z, (30) 
h(z) = K 1 {z)/z 2 + K^z)K x {z)lz + K 2 (z) - K 2 (z), (31) 



z = my/— x 2 + iO, and _Ko,i modified Bessel functions. 

Since the transverse structure is already explicit one can remove divergences in 
the right-hand side of (E9f) in an arbitrary way. For example, it is possible to continue 
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to pull out laplacians as it was down in ref. 0. However only the first term in the 
right-hand side of (|3l"|) is UV-divergent (i.e. non-integrable in the vicinity of the point 
x = 0). In fact, it is proportional to the one-loop scalar diagram so that one can also 
apply prescription (pl|) with to — 0. Note that all four terms in the right-hand side 
of (|3T| ) admit a natural diagrammatical interpretation and have simple expressions 
in momentum space, because K (z) is obtained from K 1 /z (i.e. scalar propagator) 
by the operator — m^-, and K\ just by multiplication by m 2 x 2 (i.e. — m 2 O q in 
momentum space). 

To conclude the section we note that in other theories with Abelian gauge sym- 
metry the situation is quite similar, with unessential additional considerations. For 
example, in scalar electrodynamics, one should take into account massive tadpoles 
(which are zero in QED). The tadpoles are not still covered by prescriptions (|10|)- 
(|i~2f) . However, in the next section, we shall arrive at general prescriptions, using the 
momentum space language. 

4 Differential renormalization 
in momentum space 

4.1 General renormalization prescriptions 

Let us now turn to momentum space where basic physical quantities are calculated. 
First we observe that renormalization prescriptions (p~0|) — (p~2|) are trivially transformed 
into momentum space. (We now distinguish external and internal vertices — see 
comment 5 in the end of Sect. 2.) Let Fr(q, m) = Fr(qi, . . . , q n , mi, . . . , tul) be the 
Feynman integral given by 

(2tt) 4 5 Qj) F T (q, m) = J dx x . . . dx n exp (i ^ q^x^j F T (x, m), (32) 

Fr (q, m) = J dk 1 ... dk h fl r (q,m), (33) 

where k = k±, . . . kh is a set of loop momenta of T, and Ilr is the product of propagators 
in momentum space Gi(q) = Pi(q, m) / (mf — q 2 — iO) associated with the given graph. 
Then we get the following prescriptions: 

RFp 1 = — — — — — (u/2 -D)(l- M"- 1 ) B!Fl +l ' 1 - RC y Fi +1 '\ (34) 

3 + 1 ' ' J + 1 7cr . 

A(T)Fi' l = RFp l -R'Fi' 1 , (35) 

(lo/2 - D) RF T = J2 C^RXlv = ^ RC^F T + C V F T . (36) 

7cr 7cr 

Thus the only distinction is that instead of operator S we now have dilatation 
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operator (times 1/2) 

H?4 + ^"W (37) 

By definition it acts now on integrands of Feynman integrals before integration in loop 
momenta. 

For example, for one-loop scalar Feynman integral with general masses we have 

R Idk- ln J (mf - k 2 - iO)/// 2 



(m\ — k 2 — zO)(m| — (q — k) 2 — z'O) 

dk-( — m — m 3 \ \n j+1 (mj-k 2 -i0)/fx 2 
j + 1 J 2 \ dq 1 dm\ 2 <9m 2 J (in 2 — k 2 — z'0)(m 2 — (q — k) 2 — iO) 

In particular, for mi = m 2 = and j = 0, this reduces to 

r \ i „ • w d f ^ ^ = - /dfe /ir^;" 2 ;^^, ■ (39) 



P-z'0)(-(g-£;) 2 -zO) 7 (-k 2 - i0)(-(q - k) 2 - iOf 



To calculate fl39|) it is better not to differentiate the integrand explicitly by oper- 

2^dq 

1 8 fl 2X 



ator \qjf- Rather, it is worthwhile to introduce analytic regularization [13[ 



dA J 2 H dq (-k 2 - zO) 1+A (-(g - k) 2 - z'O) 



(40) 

A=0 



When A 7^ 0, we may use the following order: to calculate the integral, differentiate 
in q, differentiate in A and finally put A = 0. When calculating the integral one uses 
the four-dimensional one-loop formula 

d * k (-k*-i0)*(-(q-k)*-i0)>° = X2 \-q 2 -i0)^- 2 > (41) 

where G is the four-dimensional G-function 

r(A 1 + A 2 -2)r(2-A 1 )r(2-A 2 ) 

G(Al ' A2) " F(A 1 )F(A 2 ) F(4-A 1 -A 2 ) " (42) 

In particular, 

G(U + A) = A(^A)- (43) 

Finally we have 

R I ^ (-f-^ -*)'-«) ^ " 2 ( : " " !0)/ " 2 ) ■ (44) 

For example 2 of Sect. 2 with a nest of divergent subgraphs, as well as for other 
examples, momentum space versions are quite similar, with product of propagators in 
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the Feynman integral in momentum space and operators D (instead of — S) associated 
with external momenta and internal masses of corresponding subgraphs. 

Note that the arguments used to prove the coordinate space prescriptions can 
be also translated in momentum space language.0 In fact, one starts with an in- 
completely renormalized Feynman integral and writes down the formula of (strictly 
four-dimensional) integration by parts. Then the coordinate space procedure of ex- 
tension of the product of distributions to the whole space looks, in momentum space, 
as dropping surface terms in this formula (which are polynomials in external momenta 
and correspond to counterterms) . To resolve the structure of (generally overlapping) 
divergences one here uses sectors in momentum space. 

Let us now remember that we considered our prescriptions (|Il3l)-(|l2|) and their 
momentum space versions (|34])-(|36|) only for graphs that do not contain tadpoles and 
that do not lead to tadpoles when inserting associated counterterms. It turns out that 
the momentum space prescriptions already have a desired form that is applicable for 
general graphs. In particular, for the tadpole graph shown in Fig. la, the recipe 





a 



Figure 1: (a) tadpole; (b) sun-set diagram 

(that is obtained using the above arguments based on integration by parts in 
momentum space) gives 

dfc /■ / 1 9\. , An(m 2 - k 2 - zO) / fi 2 

R / — = = / dk 1 - -m— ) (1 — ilfi 



m 2 — k 2 — iO J \ 2 dm J v "' m 2 — k 2 — iO 

m2 _p_ zQ)2 {t- 2 " k2 HI + m 2 /(-k 2 - ,0))] 
-2m 2 ln(l + m 2 /(-k 2 - z0)) + (m 4 /k 2 ) \n(-k 2 - i0)/fi 2 } . (45) 

If we spoil the tadpole by ln(— k 2 — iO) / /i 2 , rather than by ln(m 2 — k 2 — iO) / [i 2 we get 
a more simple result: 

r dk = r dkln(-k 2 -i0)/fx 2 
J m 2 -k 2 -i0 J {-k 2 -i0)(m 2 -k 2 -i0) 2 ' { } 



4 This was done in ref. fll^ ] where more general prescriptions were formulated for logarithmically 
divergent diagrams with simple topology of subdivergences. 
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Let us now consider the sun-set diagram shown in Fig. lb. If we used coordinate 
space arguments and started from the incompletely renormalized diagram 



R'U r = (1 + A(12) + A(23) + A(31))n r 

(where A(12),. ..are counterterms associated with three overlapping one-loop sub- 
graphs) and wanted to extend this functional from the space with deleted origin, 
x — 0, to the whole space, we would observe that the one- loop counterterms still 
vanish once we have a; / 0. Thus we do not have enough space to perform renor- 
malization of the sun-set diagram (and other similar diagrams) in two steps because 
extension to all x requires simultaneous introduction of the overall counterterm as 
well as counterterms for subgraphs. 

To overcome this complication and arrive at general prescriptions we used, in |J, 
a trick of ref. |TJ| based on Fourier transform with respect to the squares of masses 



which were treated as squares of two-dimensional vectors. A better solution of this 
problem is just to use momentum space prescriptions fl3~4])-(|3T)D that happen to be 
more flexible. We state that these very prescriptions (|34])-(|3~6|) are valid for arbitrary 
diagrams. For example, in the case of the sun-set diagram, we have (for defmiteness, 
we have chosen the first line to 'spoil by a logarithm' for renormalization of the whole 
graph) 

RF r = (l- D q>m )(l - MlJR'F*' 1 - c 23 R J dk Hm ^~_ ^'^^ ■ (47) 

Here F^' 1 differs from by the additional factor hi(m\ — pi — iO) / in the Feynman 
integral (pi is the momentum of this line), the renormalized value for the tadpole 
(with a logarithm) corresponding to the reduced graph T/{23} is given by 

r l ik ^ m \- k l~ ^ = 1 fa fl m #- _ ,) (1 _ M . > 2 K - f 



m 2 -k 2 -i0 2 J \2 dm J" m ' m 2 - k 2 - iO 

(48) 

with m = mi, and the constant C23 is found from eq. (|36|), i.e. in our case it is 
proportional to DRF23 (in fact it does not depend on the renormalization scheme 
- compare 0). Finally, counterterms for the given graph spoiled by the logarithm 
(that contribute to R' F^' 1 ), are found from eq. (|38|) . 

4.2 VVA anomaly as an example 

Since differential renormalization is strictly four- dimensional it looks preferable for 



application in each situation when dimensional renormalization |15] meets difficulties, 
for instance, in theories with chiral and super symmetries (see, e.g. Pfl where the 
initial version of differential renormalization |l| was successfully applied in such cases, 
in particular for calculation of anomalies). It is well-known that within dimensional 
renormalization the origin of anomalies turns out to be an inconsistency in definition 
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of 75 in dimensional regularization, and strictly in four dimensions one does not 
have such problems at all. Thus it is natural to ask where do the anomalies come 
from. As it was demonstrated in || , the anomalies within differential renormalization 
appear because a system of linear equations for finite counterterms to satisfy all Ward 
identities turns out to overdetermined. 

To see that the presented version of differential renormalization is calculationally 
simple and well suited for dealing with chiral problems let us once again consider 
the calculation of the ABJ triangle anomaly [16j as an example. Let T a/ 3 p (q,p) = 



S a p P (q,p) + Sf3 ap (—q — p,p) be the sum of two triangle diagrams contributing to 
one-loop Green function of the axial current J$ p = ip^p^ip, and let p correspond 
to the axial current and q, —q — p to other two external lines. In accordance with 



general prescriptions (namely eq. (|15|) translated into momentum space language), 



the differentially renormalized quantities are given by 

[' d/c ■■ 

RjS aPp (q, P ) =-iJ j^y[(D q , p , m - 1/2)(1 - M°„) ln(m 2 - p) - *0)//i 2 




where m? — iO prescriptions are omitted for brevity. Here j = AV, VV, V A and pj 
can be chosen as the momentum of the corresponding line (between vector or axial 
vertices). According to remark in subsect. 2.2 one can use an appropriate linear 
combination R = J2j an d we will do this below. 

Since we want to have Ward identities in both vector channels, 

q a RT af3p (q,p) = (q p + Pp)T a/3p (q, p) = 0, (50) 

we should introduce the logarithms symmetrically in AV and VA lines, namely, choose 
£av — £,va in the above linear combination. 

To calculate q a RjT a p p (q,p) we use simple commutation relations (|^) and in- 
troduce auxiliary analytic regularization (as in an example in subsect. 4.1). by 

ln(m 2 — p'j — ^0)//i 2 — > ~ ^ (m 2 -p?-io) x • Then we observe that the action of the 

operator (D q ^ m — 1/2) (1 — M° pm ) reduces to calculation of the finite part of the 
Laurent expansion in A (actually the pole part turns out to be zero): 

q a RjS a p p (q,p) 1 



A=0 



(2tt) 4 (m 2 - p) - iOy 

f 1 1 1 

Q ylp^ 5 m _ ^ a m _ j£_ m _ ^_|_ 

Now we apply momentum space version of fl27|) , 

11 _ 1 

m— i> m— 4— d m— i>— d m- 



(51) 



(52) 
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From coordinate space considerations, one sees that the result should be polyno- 
mial of the second degree in q,p,m. Furthermore, well-known formulae for traces of 
products of gamma matrices show that one can put m = 0. Then one uses relation 
tr7 M 7 ! ,7«7 ( g75 = eqs. (pi), (gD, 



/ 



(-P-*0) A i(-(g- k) 2 -t0y 



m 2 G (1) (A 1 ,A 2 



Q 



(-q 2 - i0) A i+ A 2- 



with 



G (1) (Ai,A 2 



and the value fl43|) as well as 

{G(l + A,l) -2G«(1 + A,1)} 
The final result is 



r(A! + A 2 -2)r(3- Ai)r(2- a 2 ) 
r(A 1 )r(A 2 ) r(5-Ai-A 2 ) 

1/2, {G«(2,A)} 



(53) 



(54) 



A=0 



A=0 



q a RAvT a /3 P — q a RvAT c 



a/3p 



8vr 5 



q a R vv T ( 



a/3p 



4vr 2 



1/2. (55) 



(56) 



so that we have Ward identities in the vector channels (BO) for R = Rav + Rv a~ Rvv ■ 
Within this choice of renormalization one can perform calculation of p p RS a (3 P (q, p) 
using 



1 



1 



1 



m— p 1 m— p 1 — $ 



75" 



m- 



m- 



1 1 1 

-7 5 - 2m -7s- 



m- 



m- 



(57) 



instead of (0), and the same technique as before, in particular formulae (0), (f42|), 
), (|55|) , with the result 



p p RT af3p = 2miRT 5al3 



2tt 



(5? 



where T 5a p is one-loop contribution to the Green function of the pseudovector current 
* < 07p75'0- This leads to the well-known value of the VVA anomaly. 



5 Conclusion 

The prescriptions of differential renormalization scheme presented above are applica- 
ble for arbitrary diagrams and look more simple as compared with previous versions. 
Let us stress that the main features of differential renormalization are pulling out 
certain differential operators and introducing a logarithmic dependence into the dia- 
grams involved. One can define it either in coordinate or momentum space (although 
momentum space turns out to be more flexible in some situations). 
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Since the mechanism of proving Ward identities in QED at the diagrammatical 
level is very transparent, it was possible to do this automatically, to all orders. Of 
course, the diagrammatical realization of non-Abelian gauge symmetries is rather 
non-trivial. One may certainly hope that the presented version of differential renor- 
malization can be useful in treating non-Abelian gauge symmetries combined with 
super and chiral symmetries strictly in four dimensions at least in lower orders of 
perturbation theory. Since the differential operators that are pulled out during renor- 
malization have very simple commutation relations with multiplication by momenta, 
the problem here is to control the logarithms that are generated by renormalization. 

Finally we note that the presented version of differential renormalization is nat- 
urally supplied with strictly four-dimensional methods of calculation of Feynman 
integrals, in particular integration by parts (I7j which is a four- dimensional analogue 



of the method of integration by parts within dimensional regularization and is 
itself based on the differential renormalization and its infrared analogue. 

Acknowledgments. I am grateful to K.G. Chetyrkin and D.Z. Freedman for valu- 
able discussions. 
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